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We report ab initio calculations of the electronic band structure, the corresponding optical spectra, 
and the phonon dispersion relations of trigonal a-HgS (cinnabar). The calculated dielectric func- 
tions are compared with unpublished optical measurements by Zallen and coworkers. The phonon 
dispersion relations are used to calculate the temperature and isotopic mass dependence of the spe- 
cific heat which has been compared with experimental data obtained on samples with the natural 
isotope abundances of the elements Hg and S (natural minerals and vapor phase grown samples) 
and on samples prepared from isotope enriched elements by vapor phase transport. Comparison of 
the calculated vibrational frequencies with Raman and ir data is also presented. Contrary to the 
case of cubic /3-HgS (metacinnabar), the spin-orbit splitting of the top valence bands at the F-point 
of the Brillouin zone (Ao ) is positive, because of a smaller admixture of 5d core electrons of Hg. 
Calculations of the lattice parameters, and the pressure dependence of Ao and the corresponding 
direct gap Eo ~ 2eV are also presented. The lowest absorption edge is confirmed to be indirect. 

PACS numbers: 63.20.-e, 63.20.dk, 63.20.D-, 68.35.bg, 65.40.Ba, 71.55.Gs, 71.70.Ej 



I. INTRODUCTION 

Cinnabar (a-HgS) is the main ore for the produc- 
tion of mercury and, in powdered form, constitutes the 
red pigment vermillion which was already used in pre- 
Columbian Peru as early as 500 BC (Chavin Empire). 
Large scale mining of cinnabar is known to have taken 
place after the conquest of the Inca Empire (1532 AD) in 
connection with the extraction of silver from low grade 
ores. It probably led to the first pre-industrial source 
of Hg environmental pollution^ Evidence for the use of 
cinnabar as a pigment is also found in Mesoamerica, dat- 
ing back to the Olmec culture (^800 BC), where it was 
utilized in ceremonial burials and for coloring beautiful 
ceramic figurines^ It was mined in Queretaro, where 
nowadays the Institute of one of the present authors 
(AHR) is located. 

In the Eastern World, China, today the main producer 
of mercury, was early using cinnabar as a pigment. The 
best known use is found in the lacquerware of the Song 
Dynasty (1000 AD)j£ It is applied, up to now, in tradi- 
tional Chinese medicine (as Zhu Sha) to treat a variety 
of ailments including colds, insomnia, restlessness and, 
less dangerously externally, skin disorders. The inter- 
ested reader will find in the web abundant therapeutic 
as well as ordering information. 

The extraction of mercury from cinnabar is docu- 



mented in Tcophrastus of Eresus' "Book on Stones" 
(^315 BC)^: "Native cinnabar was rubbed with vine- 
gar in a copper mortar with a copper pestle" , thus de- 
scribing what is probably the first mechano-chemical re- 
action. Pliny the Elder (23-79 AD), in his natural his- 
tory^, describes not only the mechano-chemical but also 
the distillation method which seems to have originated 
from Dioscorides (40-90 AD). 

Because of the existence of large cinnabar deposits, ver- 
million was also used to illuminate ancient manuscripts 
in Spain and Italy. The largest world deposits are prob- 
ably still those of Almaden in Spain but the mines were 
closed in 2002 because of the drop in the price of mercury. 
The largest producers are now China and Kyrgyztan. 

Beautiful red single crystals of a-HgS can be pried out 
of mineral samples. It is also possible to grow them from 
the elements using conventional vapor phase techniques. 
A method to grow thin films on substrates for possible 
optoelectronic applications was disclosed in 1974<£ 

Although cinnabar is the stable form of HgS under 
NTP, a zincblende-type modification, metacinnabar (/3- 
HgS), is also found in nature and can be grown in the 
laboratory by adding traces of iron (~1%) to the ele- 
ments used in the growth. The electronic and vibronic 
properties of metacinnabar have been discussed recently^ 

The present article follows the structure of Ref. 0, with 
emphasis on the differences between the crystal struc- 



ture of cinnabar and that of metacinnabar. Cinnabar, 
e.g., has three formula units per unit cell 3 - 3 - whereas 
metacinnabar only has one. This fact is responsible for 
a much larger number of ir and Raman active phonons 
in cinnabar than in metacinnabar. Their frequencies are 
investigated here. Cinnabar has two chiral (enantiomor- 
phic) modifications (space groups no. 152 (D|) and no. 
154 (Dg), primitive cell composed of two coaxial helices, 
one with three S atoms, the other with three Hg atoms). 
These modifications rotate the plane of polarization of 
light propagating along the c-axis in opposite directions 
(optical activity). They also give rise to linear k terms 
in the dispersion relations of phonons and electrons, the 
latter even in the absence of spin-orbit (SO) coupling 
(linear terms appear in /3-HgS only when SO is taken 
into account). In the present work we have measured the 
specific heat of natural and synthetic crystals of a-HgS 
with various isotopic compositions and compared the re- 
sults with ab initio calculations based on the electronic 
band structures mentioned above. Because of the uniax- 
ial crystal structure, there are two different spectra of the 
dielectric functions, depending on whether the electric 
field is polarized parallel or perpendicular to the c-axis. 
We have calculated the corresponding electronic spectra 
and compared them with unpublished experimental data 
by Zallen et al.& 

This work is organized as follows: 

Section |TT] discusses the details of the theoretical cal- 
culations, which use the VASR2 and the ABINIT set of 
codes* 1 ^ Section [IIII provides necessary experimental de- 
tails. Section IIVI presents the calculated ab initio band 
structures (with emphasis on SO and linear k effects) and 
the lattice parameters (including the bulk modulus Bq ) 
obtained by optimization of the total enthalpy. 

Section [V] displays the electronic dielectric func- 
tions obtained from the ABINIT-calculated band struc- 
tures and compares them with the results obtained by 
Kramers-Kronig analysis of the reflectivity spectra mea- 
sured by Zallen et al.& 

Section IVil is devoted to the ab-initio calculation of the 
phonon dispersion curves and a comparison of the fre- 
quencies obtained at the T-point of the BZ with Raman 
and ir measurements. Section [Vll] discusses the results of 
our heat capacity measurements and finally Section IVIIII 
is devoted to the conclusions. 



II. THEORETICAL DETAILS 

The calculations reported here have been performed 
with two different implementations of Density Functional 
Theory (DFT), VASP and ABINIT: The Vienna sim- 
ulation package VASP (see Refs. [lllil2l and references 
therein) performs ab initio electronic structure calcu- 
lations with the pseudopotential method and with the 
spin-polarized density functional theory including a self- 
consistent treatment of SO coupling. The set of plane 
waves employed extended up to a kinetic energy cutoff 



of 370 cV for a-HgS. Such a large cutoff was required 
to achieve highly converged results within the projector- 
augmented wave (PAW) schemei 13 ' 14 The semicore 5d 
and 6s electrons of Hg were included explicitly in the 
calculations. The PAW method takes into account the 
full nodal character of the all electron charge density dis- 
tribution in the core region. The exchange-correlation 
energy was initially taken in the local density approxi- 
mation (LDA) with the Ceperley Alder prescription 1 ^ of 
the exchange-correlation energy, and we also employed 
the generalized gradient approximation, GGA, with the 
PBE prescription^ We used a dense Monkhorst-Pack 
grid for Brillouin zone (BZ) integrations in order to as- 
sure highly converged results to about 1-2 meV per for- 
mula unit. We also utilized an accurate procedure in 
the calculations in order to obtain very well converged 
forces which were used for the calculation of the dynami- 
cal matrix. At each selected volume, the structures were 
fully relaxed to their equilibrium configuration through 
the calculation of the forces on atoms and the stress ten- 
sor—. In the relaxed equilibrium configuration, the forces 
are less than 0.004 cV/A and the deviation of the stress 
tensor from a diagonal hydrostatic form is less than 1 
kbar (0.1 GPa). 

In the ABINIT implementation 1 ^! 1 ^, only the va- 
lence electrons are taken into account by using the 
Hartwigsen-Goedeckcr-Huttcr rclativistic separable dual- 
space pseudopotentials2£ with an energy cutoff of 60 Ha. 
A detailed testing was performed to ensure convergence 
at the chosen cutoff. A Local Density Approximation 
(LDA) was used for the exchange-correlation Hamilto- 
nian, but some tests were also performed with the PBE- 
GGA exchange correlation Hamiltonian2I; only small dif- 
ferences were found. A grid of 6x6x4 was used to de- 
scribe electronic and vibrational properties of a-HgS in 
both symmetry groups. 

Several tests were performed with finer grids to check 
our results. After convergence was reached, an inner 
stress within the unit cell of 1 x 10~ 4 GPa remained. 

Geometrical relaxation was performed up to forces less 
than 3xl0~ 4 eV/A. The ABINIT code was used to obtain 
the optical and the vibrational properties. In the latter 
case we used the response function implementation ! 22 ' 23 
The dynamical matrices were calculated for a grid of 
6x6x3 and four different grid shifts, with a total of 83 
matrices (including the T point). After obtaining the ma- 
trices, a Fourier interpolation was used, as implemented 
in the set of ABINIT codes (anaddb) described in Ref. 
124 , to increase the density of q points. However, in order 
to calculate the optical properties we used the ABINIT 
implementation in two different versions: (a) calculation 
of the electronic interband transitions, as described in 
Ref. l25j and (b) a different approach, where an inver- 
sion of the dielectric matrix is performed as described in 

Refs. mm 



III. EXPERIMENTAL 



Natural cinnabar crystals on dolomite were purchased 
from a mineral dealer (Sebastian Barzel, Ganghofcrstr. 
1, Berlin. Origin: China). A chemical microanalysis per- 
formed on a set of selected crystals gave an Fe content 
of less than 50 at-ppm Fe and a ratio of Hg and S of 
1:0.994(4), very close to the ideal stoichiometry. Artifi- 
cial cinnabar crystals were grown from a-HgS (Chcmpur, 
99+%) by chemical vapor transport using iodine as trans- 
port agent (~1 mgcm -3 ) in a temperature gradient of 
340°C and 180°C^ This procedure lead to the simul- 
taneous deposition of the red and the black modifica- 
tion. Well shaped crystals of the red cinnabar were care- 
fully separated and selected from the growth product. 
34 S (enrichment 99.6%) was purchased from (Trace Sci- 
ences International Corp., 40 Vogcll Road, Unit 42, Rich- 
mond Hill, ON) and reacted in the gas phase at 300°C 
with elementary Hg (Chcmpur, 99,999+%) in a ratio 1:1. 
The resulting red cinnabar was subsequently refined by 
a gas phase transport reaction similar to that described 
above. Isotope enriched Hg (isotope composition 196 Hg: 
3.43%, 198 Hg: 43.0%, 199 Hg: 30.1%, 200 Hg: 8.6%, 201 Hg: 
4.39%, 202 Hg: 8.48%, 204 Hg: 1.97%) was obtained as 
HgO (also from Trace Sciences International Corp.) and 
reacted with S at 400°C. The resulting black modifica- 
tion of HgS was subsequently refined, recrystallized and 
converted into a-HgS by a gas phase transport reaction 
as described before. The heat capacities were measured 
on samples of typically ~20 mg between 2 and 350 K with 
a physical property measurement system (Quantum De- 
sign, 6325 Lusk Boulevard, San Diego, CA) as described 
in detail in Ref. [29|. Unfortunately, for the Hg isotope 
enriched sample we could only reliably separate several 
smaller crystals the total mass of which did not exceed 
~2.5 mg. The bulk of the Hg isotope enriched sample 
consisted of black material, presumably mostly /3-HgS. 



IV. ELECTRONIC BAND STRUCTURE 

Some of our band structure calculations included SO 
interaction, others not. The latter (NOSO) were per- 
formed in order to examine the effect of the SO interac- 
tion (not always small because of the large SO coupling 
in the Hg atoms) and also to reduce computational time 
whenever needed (i.e., for the calculation of optical and 
phonon spectra). 

We show in Fig. [T]the electronic band structure of a- 
HgS in an extended energy range along several high sym- 
metry directions of the BZ (the BZ and the notation for 
the high symmetry points can be seen in Ref. l30j ) includ- 
ing SO interaction, although in this figure its effects can 
only be clearly observed for the 30 bands derived from 
Hg 5d states ( -5 to -8 eV range). Beside these bands, the 
figure includes bands in which the 3s electrons of sulfur 
are dominant (-12 to -14 eV), sulfur 3p valence bands (-5 
to eV) and empty conduction bands in which Hg 6s 




FIG. 1: Electronic band structure (VASP LDA SO) of a-HgS. 
The arrows near F are discussed in Fig. [2] 



and 6p states dominate. The calculation was performed 
with the VASP code, using the LDA to approximate ex- 
change and correlation. Similar results were obtained 
with the ABINIT code and thus will not be presented 
here. Calculations, performed with the WIEN2K code 
using the generalized gradient approximation (GGA) for 
exchange and correlation have been already published^ 
They were used in Ref. |3l| to investigate the pressure 
dependence of the band structure and the phase transi- 
tion from cinnabar to the rock salt structure (at about 
25 GPa£2). Figure [1] and [2] show an indirect valence- 
to-conduction energy gap at 0.82 eV (from near A to 
the T point of the BZ) followed by a direct one at 0.89 
eV. These gaps are considerably smaller than the exper- 
imentally reported ones (~ 2.25 eV— ). The difference is 
related to the so-called gap problem which appears when 
local density functionals are used in the calculations. It 
is semiempirically removed by ad hoc adjusting the gap 
with a scissors operator^ We present in Figs. [T] and [2] 
raw data obtained using the LDA. For the calculation 
of optical spectra we shall add 1 eV as a scissors op- 
erator, a correction determined for /3-HgS by Svanei^ 
with the GW Hamiltonian. The Hg 5d bands (cf. Fig. 
[J) cluster in two groups centered between 6 and 7.8 eV 
below the top of the valence bands (TOV) , the splitting 
resulting from a combination of the orbital field and the 
SO interaction. Correspondingly, two peaks are seen in 
the photoemission spectra centered at 7.0 and 8.8 eV4I 
While the width of these bands agrees with the calcu- 
lated one, their higher binding energy (by about ~1 eV) 
may be due to an LDA effect similar to that which causes 
the gap problem^ In Ref. [38| this energy difference was 
calculated to be 1.2 eV for wurtzite InN. The sulfur 3s 
bands are centered around 12.5 eV. They probably cor- 
respond to a peak observed in photoemission at the same 
energy (see Fig. 7 of Ref. I37I ). 

In order to better see the indirect (E{ ) and direct (Eq) 
gaps, as well as the SO splitting of the TOV, we show in 
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TABLE I: Cell volume, lattice parameters, u and v and bulk moduli of a-HgS resulting from calculations compared to experi- 
mental data. 
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FIG. 2: Electronic band structure (VASP LDA SO) of a-HgS 
in an enlarged energy scale. Eq and Ei mark the direct and 
the indirect gap, respectively. 



Fig. [5] the bands of Fig. [T]in an expanded energy scale. 
A lowest indirect gap Ei = 0.82 eV is observed in this 
figure, followed by the already mentioned direct one at 
0.89 eV. The SO splitting at T is found to be |A | = 39 
mcV. An interesting question is the sign of this splitting 
in view of the fact that, anomalously, that of /3-HgS has 
been found to be negative because of admixture of Hg 5d 
electrons to the dominant S ip electrons. Depending on 
the amount of admixture, the splitting, which is rather 
small anyhow, can become either positive or negative.— 
We shall present evidence that the spin-orbit splitting An. 
of a-HgS is positive, contrary to that of /3-HgS. 

By downloading the eigenvectors obtained from the 
VASP code, with and without SO coupling, we have ob- 
tained the symmetries of the corresponding states around 
the direct gap (at T). The symmetries corresponding to 
the simple D3 point group (without spin) are two singlets 
(ri and T2 ) and a doublet T3.— The double group (in- 
cluding spin) has three representations with symmetries 
r 6 (a doublet), and two singlets (r 4 , F5), degenerate on 
account of time-reversal.— The symmetry of a single spin 
is T§. We found the state at the bottom of the conduction 
band (Fig. [5]) to have Ti symmetry, mainly consisting on 



each atom mainly of 6s and Qp z ) atomic functions of Hg 
(see Fig. [3]) . Its symmetry becomes r 6 when SO coupling 
is taken into account. The top of the valence band at F is 
an orbital doublet ( T3 symmetry), composed mainly of 
(3p x , 3p y ) functions of sulfur, with some admixture of 5d 
functions of Hg. The latter are responsible for the small 
(in /3-HgS negative) value of An. The next lowest T state 
(0.4 eV below T3) has T2 orbital symmetry. Under SO 
interaction, the T3 state splits into two spin doublets, the 
upper one having (r4, Y§) symmetry, the lower Tq . This 
ordering corresponds to a positive SO splitting, in which 
for the upper state the angular momenta of the spin and 
the orbit are parallel (see Fig. 5.1 of Ref. 1421) Hence these 
calculations suggest that An is positive for a-HgS. 
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FIG. 3: (Color online) Electronic band structure (VASP LDA 
SO) of a-HgS with partial electronic DOS as labeled in the 
inset of the r.h.s. of the figure. 

Time reversal compels all double group states to be 
doubly degenerate at T. This degeneracy, however can 
split away from T, resulting in the so-called linear T 
terms which we have already discussed for /3-HgS (zb 
structure) .— While in /3-HgS only r 8 quadruplets exhibit 
linear in k splittings, most states (not all) show these 
splittings for the a-HgS structure. Most interesting is 
the splitting of the T3 orbital doublets for k along the 



trigonal axes (T - A in Figs. [2] and [3]). This splitting 
is related to the chirality of the structure and its op- 
tical activity: it reverses its sign when going from the 
space group no. 152 (D|) to the space group no. 154 
(Df) structure, as does its optical rotatory power 4^ The 
structure of the linear in k splittings can be investigated 
using so called invariants, i.e. Hamiltonian-likc functions 
linear in k which have the full symmetry of the crystal^ 
Such invariant is, for the symmetry at hand: 



Hf. — Lj z J z k z + C X {J X k x + Jy ky), 



(1) 



where the coefficients C depend on the state under 
consideration and J transforms like angular momentum 
J = 1 under the symmetry operations of D3 . In the 
absence of spin, the expectation values of Eq. (TTJ) vanish 
for states of symmetry Ti, T2, r 4 , and T$ . As already 
surmised, the only coefficient which does not vanish in 
the absence of spin is C z . We have calculated for the 
top of the valence band in the absence of spin C z = 0.22 
eV A. When spin is included, all states of a-HgS exhibit 
splittings linear in k except those of (T4, T^) symmetry 
for k along A-M. This is shown in Fig. @] for the lowest 
conduction band and the 3 highest valence bands at T. 
The calculated values of the coefficients Ck (the splitting 
is defined as (±) Ck • k) are displayed in Table Hill Note 
that for k along T-M and T-K Ck is small but not exactly 
zero for the highest valence band. This is probably a 
computational error resulting from having had to deal 
with very small energy differences. 



0.588 — 
0.584 
0.58 — 
0.576 



>■ 

at 

<5 -0.32 — 
C 
LLI 



I 

0.001 



T 




M^r^A 



r->K 



FIG. 4: Magnification of the electronic bands showing the lin- 
ear splitting of the lowest conduction and the highest valence 
band near F (in units of 2vrxA _1 ). 

It is easy to understand this vanishing of Ck for k along 
r-M and T-K on the basis of Eq. ([!]). The T 4 and T 5 dou- 
blet has the same transformation properties as J z =3/2 
and J 2 =-3/2 angular momentum states. J x and J y can 
only connect states which differ in angular momentum 



J z by ±1, hence there can be no linear splittings of the 
T4 , r 5 states for k along either k x or k y . This observa- 
tion, represented also in Fig. |H indicates that the (F4, 
r$ ) states correspond to the top of the valence band and 
confirms the positive sign attributed to An. above. No- 
tice that the non-zero values of Table ITU especially those 
for k along the T-A direction, are similar to those found 
for /3-HgSii However, while a sign was attached to the 
values in Table [TT1 of Rcf. 0, corresponding to placing the 
anion at the origin of coordinates^, the intricacies of a- 
HgS and the corresponding VASP program did not allow 
us to make this assignment in Table [TT] with respect to 
the two possible enanthiomorphic varieties (the signs for 
space group no. 152 variety should be opposite to those 
for no. 154). This sign should determine the correspond- 
ing sign of the dispersion of the optical activity, which 
has hitherto not been determined.— 



V. OPTICAL PROPERTIES 

We display in Fig. [5] the calculated real and imagi- 
nary parts of the dielectric function for the electric field 
perpendicular and parallel to the trigonal axis for a-HgS. 
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FIG. 5: (Color online) Calculated spectra (ABINIT LDA, 
NOSO) of ei,2 (□ and o, respectively) for electrical field (a) 
E ± c and (b) E \\ c. 

The calculations were performed using the ABINIT 
code without SO coupling and with the LDA Hamilto- 



A dA /dP dlnAo/dlnV E dE /dP dEo/dlnV 
(meV) (meV/GPa) (meV) (meV/GPa) (meV) 



Q-HgS VASP-LDA 39 

/3-HgS VASP-LDA -111 
a-HgS ABINIT-LDA 

/3-HgS ABINIT-LDA -180 
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TABLE II: Spin-orbit splitting Ao at the top of the valence band of cinnabar (a-HgS) and metacinnabar (/3-HgS) and its 
derivatives vs. pressure and volume (the latter logarithmic). Also, direct gap at the P-point of the BZ (not corrected for the 
"gap problem") together with the corresponding pressure and volume derivatives. 



energy C (T - A) C (T - M) C (P - K) 
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TABLE III: Coefficients of the terms linear in k for the states 
around the direct gap of cinnabar. The linear splitting is ±C. 
These splittings were calculated including spin-orbit interac- 
tion. Because of the chirality, linear terms in k also appear in 
the absence of spin-orbit interaction for orbital states of T3 
symmetry such as those at the top of the valence band. In 
this case C = 1.357/2tt = 0.216 eVxA. 
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man for exchange and correlation. We separated the cal- 
culated conduction and valence bands by a leV "scissors 
operator" which we took to be equal to that obtained 
with the GW approximation for /3-HgS^ A Lorentzian 
broadening of 0.1 cV was applied to all transition en- 
ergies. Notice that ei becomes zero for E ps17 eV, an 
energy which corresponds to the plasma frequency of the 
valence electrons, with some contribution of the 6d semi- 
core electrons of Hg. This plasma frequency can be ob- 
tained more accurately from a plot vs. energy of the 
electron energy loss (EEL) function 



F(EEL) = -9f(l/e), 



(2) 



as shown in Fig. [6] Both EEL functions, for E paral- 
lel and perpendicular to the c-axis, show a broad band 
peaking around 17.5 eV. 

The analysis of the overall optical response in the 0-20 
cV range can be completed by calculating the effective 
number of electrons contributing up to the energy E, 
using the expressior 



,46 



N, 



off 



&£«**'* 



■k 4ire' 



(3) 



The energy dependence of N e f[ obtained by substitut- 
ing into Eq. ((3]) the values of e\{E) given in Fig. [5] are 
shown in Fig. [7] According to this figure, N e g has not 
attained saturation at the maximum energy of «19 cV. 



FIG. 6: (Color online) Calculated electron energy loss func- 
tion function of a-HgS. 



The values reached at this energy are 5 electrons per 
atom for the parallel component and 4 for the perpen- 
dicular, rather close to the average number of valence s 
and p electrons per atom. Hence we conclude that the 5d 
electrons of mercury contribute to N e g mainly at much 
higher energies^ It is interesting to note that up to 19 
eV N e s is larger for the parallel (5 electrons/atom) than 
for the perpendicular (4 eV/atom). It is not easy to find 
a simple reason for this fact, which is also reflected in 
the higher value of £i(0) for E\\c (ordinary rays) than for 
E _L c (extraordinary rays) see Fig. [S] 

The calculated spectra of Fig. [S] arc compared in Fig 
[5] with unpublished experimental results of Zallcn et al. 
obtained by Kramers-Kronig analysis of normal incidence 
reflectance data£ The experimental values are through- 
out the whole spectral range lower than the calculated 
ones, a fact which cannot be completely attributed to ex- 
perimental errors: measurements at energies below the 
gap (E ps 0) give t\ = 8.4 for E || c and ei = 6.9 for 
E _L c (note that ei = 8.4 is the highest value of ei(0) 
found for any mineral). However, although the calcula- 
tions give values of e larger than the experimental ones, 
the calculated birefringence (exii - ei_i_) and its sign agrees 
roughly with the measured one4£ The calculated ei and 
62 show considerable structure: we have indicated the 
most prominent peaks in 62 by arrows. Most of this 
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FIG. 7: Effective number of electrons per atom N e s contribut- 
ing to the optical absorption up to the energy E, as obtained 
by integrating the e^ of Fig. [5] with Eq. JT}. 



FIG. 8: (Color online) Calculated (ABINIT LDA, NOSO) (□ 
and o) and measured spectra (black and red solid lines) of 
61,2 (E _l_ c) and ei^ (E || c). The experimental data have 
been taken from Zallen et at. (Ref. @). The vertical arrows 
mark structures in the calculated e(w) curves. Their energies 
are listed in Table ITVl 



structure seems to correspond to features observed ex- 
perimentally: we have listed in Table IIVI the calculated 
positions of peaks in ei and compared them with exper- 
imental ones. The calculated absorption starts at about 
2 eV but this energy depends on the value of 1 eV, some- 
what arbitrarily taken for the scissors operator. Accord- 
ing to Zallen 3 -^ an indirect gap appears about 10 meV 
below the direct one. In Fig [2] such gap can indeed be 
seen about 40 meV below the direct one. The difference is 
not surprising as it may be related to the gap problem al- 
ready mentioned. In any case, the direct-indirect energy 
differences are so small that even ab initio calculations 
are not likely to reproduce them exactly. It is already 
remarkable that we can get its sign. Notice in Fig. [2] 
of Ref. [3a that near the absorption edge the extraordi- 
nary rays (E _L c) are more strongly absorbed than the 
ordinary ones (E \\ c). This is the opposite of what is 
shown in Fig. [5] for the overall absorption spectrum and 
its main peak at about 4.5 eV. Hence a crossover of the 
corresponding values of 62 must take place at energies 
slightly above the absorption edge. Unfortunately in this 
region the calculations, as well as the extant measure- 
ments, are not accurate enough to reveal this fact. 



VI. PHONON PROPERTIES 

Regrettably there arc no experimental data for the 
phonon dispersion relations throughout the whole Bril- 
louin zone (BZ). Inelastic neutron scattering measure- 
ments are, in principle, possible. However, the rather 
large absorption cross section for the isotope 199 Hg (~ 
17% abundance in the natural isotope mixture) would 
require an isotope enriched sample. Given this problem, 
the most appropriate technique to map the whole BZ 
would be inelastic x-rays scattering. Although the resolu- 
tion of this technique is rather limited (~10 cm -1 ), small 
crystals ( ~ 1 mm 3 ) can be handled. At this point, the 
only available experimental data are those obtained with 
Raman and ir spectroscopy ! 48 ! 49 In Ref. |49| semiempiri- 
cal calculations are presented based on a mixed valence- 
Coulomb force field adjusted so as to fit as well as pos- 
sible the experimental results of Ref. |48|, consisting of 
Raman and ir- active phonons at or very near the cen- 
ter (r) of the Brillouin zone. The 6 atoms per primitive 
cell give rise to 18 vibrational modes, 3 of which have 
zero frequency at T (acoustic modes). Thus we are left 
with 15 modes, 5 T 3 doublets (ir and Raman active), 
2 Ti singlets (Raman active) and 3 T2 singlets (ir ac- 
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4.94 
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5.40 


6.47 7.17 
6.61 - 



TABLE IV: Comparison of calculated (see Figs. [8ja) and (b)) and experimental energies (in eV) of peaks in £2(1^) for E _L c 
and E II c£ 



tive)4£ The ir active modes split into longitudinal and 
transverse, depending on whether the E'-ficld is parallel 
or perpendicular to the scattering vector, the splitting 
being determined by Born effective charges4£ 
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339 
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a Ref. 49 



TABLE V: Phonon frequencies measured for a-HgS. Also, 
semiempirical valence-Coulomb force field calculations by 
Nusimovici and Gorre.— and the results of our ab initio calcu- 
lations with spin-orbit coupling (SO) and without spin-orbit- 
coupling (NOSO) are given for comparison. 

The ratios of the L to T frequencies of the ir-active 
modes (1^ and T3) can be used to calculate the corre- 
sponding ratio of eoAoo dielectric constants using the 
generalized Lyddane- Sachs- Teller (LST) relation: 



We performed ab initio calculations of the phonon fre- 
quencies vs. wavevector k using the ABINIT code within 
the LDA, as described in Section |TTJ We had realized in 
past work that the calculated phonon frequencies could 
depend considerably on SO interaction in materials con- 
taining heavy ions (Bi, PbTe, HgTe)<^r— Since mercury 
falls in this category, we performed a few calculations 
of phonon frequencies taking SO interaction and the re- 
laxed lattice parameters (cf. Table [TJ) into account. Ta- 
ble [Vj summarizes the results of such calculations for the 
phonon frequencies at the T point, together with their 
counterparts obtained without SO interaction (NOSO). 
For comparison, experimental data from Raman and ir 
spectroscopy^ and the results of semiempirical calcula- 
tions fitted to the experimental data are also listed4^ The 
effect of SO interaction on the calculated frequencies is 
small, typically 1% less than when it is not taken into ac- 
count. Most of the calculated frequencies become slightly 
closer to the measured ones than in NOSO calculation. 
Since the discrepancy between calculated and measured 
frequencies lies between 5 and 10%, we decided, however, 
not to include SO interaction any further in the ab ini- 
tio calculations, so as to reduce computational time. The 
calculated LO-TO splittings of most modes are somewhat 
smaller than the measured ones, a fact that had already 
been observed in earlier work on zincblende and rocksalt 
type materials i£i On the whole, the agreement between 
measurements and calculations is reasonable, especially 
when one considers the complicated crystal structure of 
a-HgS, with six atoms per unit cell. 

We display in Fig. [9] the phonon dispersion relations 
calculated with the NOSO restriction. It is worth notic- 
ing in this figure that the five T3 doublets at the T point 
split linearly in k along the T-A direction, an effect sim- 
ilar to that discussed in Section [TV] in the absence of SO 
interaction. These splittings, induced by the chirality of 
the structure, should be responsible for optical activity 
related to the ir-active plionons. 



n^y = -- 

U? too 



(4) 



From the calculated NOSO frequencies in Table [V] we 
obtain using this LST relation: 

(eo/eoo)||= 1-20 and (e /eoo)± =1-24. Using £ m from 
Fig. [5] and the values of eo from the ABINIT pro- 
gram (e || = 12.39 and eoj_=10.33) we find the ratios 
( e o/ e oo)|| = l-23 and (eo/ e oo)± =1.23 in good agreement 
with those obtained with the generalized LST relation. 



We display in Fig. [TOTa) the phonon density of states 
(DOS) calculated from the dispersion relations of Fig. 
[HI together with its decomposition in its S-like and Hg- 
like partial components. As expected, the low frequency 
band (0 - 100 cm -1 ) corresponds mainly to Hg vibra- 
tions whereas that between 230 and 350 cm -1 is mainly 
sulfur-like. It is interesting to notice that the interme- 
diate band, between 110 and 190 cm -1 , is almost pure 




FIG. 9: Calculated (NOSO) phonon Dispersion of a-HgS. 



sulfur-like. We have plotted in Fig. ITUTb') the optical two- 
phonon DOS calculated with the restriction that the sum 
of the two participating phonon wave vectors be zero. 
Both DOS, corresponding to sums and differences of two 
phonons are given. These DOS will be of interest in the 
interpretation of carefully taken ir and Raman second 
order spectra, a rather delicate task in view of the large 
number of Raman and ir phonons active in first order. 
The difference spectra should vanish at low temperature 
according to the well known Bosc-Einstein statistical fac- 
tors. This, as well as the corresponding temperature de- 
pendence of the sum spectra, should help in identifying 
the second order Raman structures. 




100 200 300 400 

E (cm 1 ) 



o 






sum 

difference 




(b) 


|1.5 








(states/ 

b 








CO 

O 0.5 

Q 










.J i .Li i 


k 



200 400 600 800 

E (cm" 1 ) 



VII. SPECIFIC HEAT 

Heat capacity measurements of a-HgS have been per- 
formed before by Khattak and coworkers in the temper- 
ature range up to ^70K^ In order to compare with our 
theoretical calculations we have extended these measure- 
ments up to ^350K on a set of samples comprising nat- 
ural and artificial crystals of ev-HgS. The results of these 
measurements are displayed over the whole temperature 
range in Fig. [TT] 

To emphasize the low temperature regime our experi- 
mental data are also shown in the standard C p /T 3 plot 
('Dcbyc plot') together with the data obtained by Khat- 
tak el al. for comparison (Fig. IT^|) . Additionally, the 
results of our calculations based on the standard inte- 
gration of the DOS of Fig. [TUTa) with the appropriate 
Bose-Einstein factors are also provided in both figures. 
Fig. [T2] illustrates the very good agreement of our exper- 
imental data with those of Khattak el al.. Furthermore, 
there is also no noticeable difference between our natural 
and artificial a-HgS sample. 

Figure [T^] reveals the characteristic maximum in the 
quantity C p /T 3 at ~10K which had actually been ob- 
served at 7K for the zincblcnde modification (/3-HgS). 
The position of this maximum is usually determined by 
the lowest maximum of the phonon DOS which, accord- 
ing to Fig. ITUTa) takes place at 65K. The ratio 65/10=6.5 



FIG. 10: (Color online) Phonon DOS including partial con- 
tributions, Hg and S, also two-phonon DOS for sums and 
differences of two phonons with equal k-vectors are shown. 




o a-HgS 

LDA-ABINIT NOSO 

LDA-ABINIT NOSO + 1.88xlo J J moF 1 K 2 xT 



200 

7(K) 



400 



FIG. 11: (Color online) Heat Capacity of the artificial sample 
of a-HgS (o). The results of our calculations (C V (T), NOSO) 
are represented by the (red) solid line. The (black) solid line 
shows C P (T) obtained by adding a term linear in temperature 
as described in detail in the text. The (black) dotted vertical 
marks the Petit-Dulong limit of 6 R, where R is the molar gas 
constant. 
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is typical for the maximum of C p /T 3 in many semicon- 
ductors. Our theoretical data agree fairly well with the 
experimental data. Near the maximum in Cp/T 3 the- 
ory falls short by about 3.5%. Above 20 K, up to 250 
K there is good agreement of experimental and theoret- 
ical results. Above 250 K the theoretical data approach 
the Petit-Dulong value of 6 x R, where R is the molar gas 
constant, somewhat slower than the experimental data. 
This difference is due to increasing importance of the 
difference of the constant-volume and constant-pressure 
specific heats which is given by^: 



where 



C P (T) = C V (T)+A(T)-T, 



A{T)=a 2 v {T)-B -V moh 



(5) 



(6) 



a v (T) is the temperature dependent coefficient of the vol- 
ume thermal expansion, Bq the bulk modulus and Knoi 
the molar volume at T = K which amounts to 28.5 
cm 3 mol -1 . 

Good agreement of the experimental and theoretical 
data is achieved if we choose for the coefficient of the 
linear temperature increase an average value Aq of 

A(T) « A = 1.88 x Kr 5 Jmol _1 K -2 . 

For an uniaxial system the volume coefficient of ther- 
mal expansion can be calculated from the linear thermal 
expansion coefficients, a» and a±, according to 



a v (T) « ai|(293J0 + 2 • a ± {293K). 



(7) 



Using for the linear expansion coefficients the literature 
values for temperatures between 20°C and 200 °C, a\\ = 
18.8 xlO" 6 K" 1 and a± = 18.1 xlO" 6 K" 1 given by 
Ohmiyaa*, the experimentally determined linear increase 
given above implies a bulk modulus, B , of 
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FIG. 12: (Color online) C p /T 3 vs. T representation of our 
heat capacity data. Natural a-HgS sample: a; artificial vapor 
phase grown sample: o. The solid line represents our theoret- 
ical data (NOSO, for details see text). Earlier heat capacity 
data by Khattak et a/.— are represented by: o. 
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FIG 13: (Color online) Cp/T 3 vs. T representation of the heat 
capacities of a-HgS samples grown from isotope enriched Hg 
and S as indicated in the inset. 



Bo =41(l)GPa 
in good agreement with the theoretical results (cf. Table 

Like in some of our previous wor k 51 ' 52 we have investi- 
gated the dependence of C p /T 3 on the isotopic masses of 
the two components of cinnabar and compared our exper- 
imental data with the theoretical results. Fig. [TBIdisplavs 
the low-temperature heat capacities of a- 1989 Hg nat S and 
a- na *Hg 34 S. As expected, the heat capacity is slightly 
larger for the Hg isotopically enriched sample. Small dif- 
ferences in the heat capacity of the S isotope sample to 
the sample made from components with the natural iso- 
tope composition appear at higher temperature and are 
barely seen in the C p /T 3 representation. They are, how- 
ever, clearly revealed in the derivatives of Cp/T 3, versus 
to the isotope masses displayed in Fig. Q3] 



Two pronounced features are observed in the logarith- 
mic derivative versus the mass of the S constituent (Fig. 
Q3Ja)), a rather pronounced sharp peak at ^6 K and a 
broad hump with maximum at ^100 K and probably a 
shoulder at ^50 K. This broad feature reflects the two 
broad S-like bands in the phonon DOS (cf. Fig. [TU] be- 
tween 150 and 180 cm' 1 and between 250 and 350 cm -1 . 
The origin of the sharp low-temperature peak at 10 K is 
not immediately obvious from the phonon DOS. We at- 
tribute it to a small mixture of Hg- and S-likc phonons 
which also explains the small spike in the S-likc phonon 
partial DOS at about 50 cm -1 . All features, especially 
the broad hump peaking at ^100 K, are fairly well re- 
produced in position and magnitude by the results of our 
calculations. 

The agreement between experiment (two separate runs 
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on two separately prepared samples from the same batch 
of isotopically enriched Hg) and theory for the logarith- 
mic derivative w.r.t. the mass of the Hg atoms (for the 
mass of the Hg constituent we took the mass-average of 
the Hg isotopes weighted by their abundance). A peak 
occurs in the logarithmic derivative at ^6 K which corre- 
sponds well with the position expected from theory. But 
the experimental peak overshoots the theoretical calcu- 
lation by almost a factor of three. The origin of this dif- 
ference might be due to a small traces of metacinnabar 
/3-HgS or other related effects. /3-HgS has a significantly 
higher low-temperature heat capacity than a-HgS, in a 
Cp/T 3 representation the peak is by a factor of three 
higher than that of a-HgS.— Comparing the heat capac- 
ities of the two modifications of HgS we estimate that 
a 1.5% admixture of /3-HgS into our sample, probably 
hardly noticeable by visual inspection, would be suffi- 
cient to explain the increase of the magnitude of the 6 K 
peak in the logarithmic derivative versus the Hg mass of 
a-HgS. The high cost of mercury isotopes has prevented 
a more detailed analysis of this anomaly. 

Finally, we demonstrate the relationship of the log- 
arithmic derivatives versus temperature and versus the 
masses of the two constituents following several previous 
results Jii 

For the relation of the temperature and mass deriva- 
tives in a two component system we have used the ex- 
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(8) 



where Mi and M 2 are the masses of the two con- 
stituents, i.e. Hg and S, respectively, for a-HgS. 

This relation is nicely obeyed for a-HgS, as is is demon- 
strated in Fig. \T5\ Especially, the low temperature peak 
which emerges essentially from the logarithmic derivative 
versus the mass of Hg is very well reproduced, thus prov- 
ing that the differences between experiment and theory 
discussed above originate from the experimental data. 



VIII. CONCLUSIONS 

Ab initio electronic band structure techniques, espe- 
cially those which use available computer codes such as 
VASP and ABINIT, are powerful methods to investigate 
electronic, optical, vibronic, and thermodynamic proper- 
ties of crystals and compare the results with experimen- 
tal data. Here we apply these techniques to cinnabar 
(a-HgS) which has a chiral structure (space groups D\ 
and -Dg, three molecules per primitive cell) more com- 
plicated than those usually dealt with. We calculate the 
electronic band structure and the spectral dependence of 
the dielectric function, confirming the indirect nature of 
the lowest gap closely followed by a direct gap. We inves- 
tigate the linear terms in k , some of which appear even 
in the absence of spin as a result of the chirality. We 



FIG. 14: (Color online) Logarithmic derivatives with re- 
spect to the masses m of the constituents Hg or S, 
dln(C p /T' ! )/dlnm of the experimental heat capacities, (a) 
Compares a- nat Hg nat S with a sample of a- nat Hg 34 S. (b) Log- 
arithmic derivative calculated by comparing samples of a- 

nat Hg natg and Q _198.9 Hg natg ag described m detail in the text . 

Different symbols (o, □) indicate independent runs on sam- 
ples of the same preparation. The (red) solid lines represent 
the logarithmic derivatives obtained from our theoretical data 
calculated for various isotope masses. 



compare the calculated dielectric function with unpub- 
lished data by Zallen et al. . We use the ABINIT code to 
calculate the frequencies of Raman and ir phonons and 
their dispersion relations. The densities of states of one 
and two phonons are also calculated. We devote the last 
section to present experimental data on the specific heat 
vs. temperature for natural samples and especially grown 
isotopically modified ones. These results are compared 
with ab initio calculations. Generally, good agreement 
between experiment and ab initio results is observed. 
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FIG. 15: Illustration of the relationship of the logarithmic 
derivatives w.r.t to temperature and the respective isotope 
masses given in Eq. [5] The (blue) dashed and the (red) 
dashed-dotted lines represent the logarithmic derivatives ver- 
sus the isotope masses obtained from our theoretical data for 
Hg and S, respectively. The (black) solid line is the sum of 
the two mass derivatives, o denote the logarithmic derivative 
versus temperature of a- nat Hg na S. 
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